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Abstract 

We study the weak productions of novel heavy mesons, such as n' c , h c , h' c , x' c oj A"(3940), F(3940), X(3872), and 
y(4260), in the semi-leptonic B c decays. Since there is still no definite answer for the components of X(3940), 
y(3940), X(3872), y(4260) so far, we will assign them as excited charmonium states with the possible quantum 
numbers constrained by the current experiments. As for the weak transition form factors, we calculate them in 
the framework of light-cone QCD sum rules approach, which is proved to be a powerful tool to deal with the non- 
perturbative hadronic matrix clement. Our results indicate that different interpretations of X(3940) can result in 
remarkable discrepancy of the production rate in the B c decays, which would help to clarify the inner structure 
of the X(3940) with the forthcoming LHC-b experiments. Besides, the predicted large weak production rates of 
X(3872) and F(3940) in B c decays and the small semi-leptonic decay rate for B c — * F(4260) all depend on their 
quantum number J PC assignments. Moreover, the S — D mixing of various vector charmonium states in the weak 
decay of B c is also discussed in this work. The future experimental measurements of these decays will test the 
inner structures of these particles, according to our predictions here. 

PACS numbers: 14.40. Gx, 13.20. Gd, 11.55.Hx 
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I. INTRODUCTION 



A number of new hidden charm states were observed recently by experiments, such as X(3872), 
X(3940), Y(3940), Z(3930) and Y(4260) Q, Q, H, Q, H] - Their quark structures are still not fully under- 
stood 6J|. In particular, the X(3872), which exhibits various impenetrable aspects, is labeled as the poster 
boy of the new heavy hadrons [?|]. Although the quantum numbers J PC = 1 ++ of X(3872) are strongly 
favored by the experiments, there is not a definite answer on its components yet due to the fact that 
none of the interpretations can fit all the available experiments satisfactorily. The assignment of X(3872) 
as a 2 3 Pi charmonium state, even without the mass gap problem 1 as claimed by calculations based on 
the Lattice QCD recently in also bears other difficulties. The tiny decay width of AT(3872), whose 
upper bound is 2.3 MeV with 90% confidence level, is much less than the number predicted in theory 
\t\. Another puzzle is the G parity violation indicated by the measurement of the ratio of branching 
fractions ~^w^^jM^t^\ " = l-0db0.4±0.3 [9, 10, 11]. The difficulties of the charmonium interpretation 
invoke various models for the structure of X (3872), such as multi-quark state [3]) hybrid meson fl^ . 

191 ] and so on. In one word, the inner structure of X(3872) is 



15, 



16 



17 



18. 



nuclear-like molecular state 
still not settled down. 

In addition to the intriguing particle X(3872), other heavy hidden charm mesons X/Y(3940), Z(3930) 
and Y(4260) mentioned above also attract comprehensive attention recently Q], among which Z(3930) 
can be well established as the first radial excited states of tensor charmonium Xc2 reasonably and will 
be left out in this paper. Even though the experimental results of h c and r]' c are essentially consistent 
with theoretical expectations, there are still some particular aspects deserving further investigations 
20, 22, HQ. Besides, we also predict the production rate of h' c state in the weak B c decays, which 
has not discovered. In addition, the 25 — ID mixing of '0(3686) and ^(3770), which is of great interest 
in quarkonium physics, is considered in the weak decays of B c . More important, we also investigate the 
production of Y(4260) and ^(4415) in the weak B c decays as the mixing of AS and 3D states. For the 
completeness, the 35" — 2D mixing of ^>(4040) and ^(4160) in the B c decay is also included. 

In this work, we do not attempt to discuss all the explanations of these states. Instead, we concentrate 
on the assignment of these heavy mesons as charmonium states with the possible quantum numbers 
constrained by the available experiments and then study their production properties in the B c decays. 
To be more specific, we will assign the X(3872) as a 2 3 Pi charmonium, the Y (4260) as a 4 3 Si charmonium, 
the X(3940) being either 3 1 Sq or 2 3 Pi charmonium, and 2 3 Pq(cc) state for the Y(3940), the quantum 



1 The mass of 2 3 Pi charmonium predicted by the quark model is about 100 MeV larger than the measured X(3872). 
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numbers of which as the charmonium states are most favored by the current experiments [2CJ,|21[ , although 
there is overpopulation of 2 3 P\ (cc) meson in the charmonium famil y. The non-leptonic weak decays 
B + — > X(3872)K + and B c — > X(3872)tt(K) have been studied in [251 . |26I ] in order to pry information of 
the inner structure of X(3872). It is found there that different J assignment of X(3872) will give quite 
different decay rate for future experiments to measure. Here we will focus on the semi-leptonic weak 
production of charmonium particles in B c decays, where it is theoretically easier compared with that of 
non-leptonic decays. We will see that the different assignments of J PC quantum number to X(3943) 
will give remarkable different branching ratios of semi-leptonic decays. Therefore our predictions can be 
used by future experiments to test the quark structures of these mesons. The main job of calculating 
the branching fractions of the semi-leptonic decays of B c is to properly evaluate the hadronic matrix 
elements for B c — > M C c {M= pseudo-scalar (P), scalar (S), vector (V) or axial vector (A) charmonium), 
namely the transition form factors. 

The precise calculations of form factors are very complicated due to the non-perturbative QCD effects 
in the hadron as a bound state. Several methods have been developed to deal with this problem on 



the market so far, such as simple quark model 
(SVZ) 



31 



321 ]. light-cone QCD sum rules 



33 



34, 



27 



light-front approach 



28, y, 0, 



QCD sum rules 



35f | . perturbative QCD factorization approach 36l. 137]. 



Although the QCD sum rules approach has made a big success, short distance expansion fails in non- 
perturbative condensate when applying the three-point sum rules to the computations of form factors in 
the large momentum transfer or large mass limit of heavy meson decays. The light-cone QCD sum rules, 
as a marriage of QCD sum rules techniques and the theory of hard exclusive processes, were developed in 



an attempt to overcome the difficulties 

y,y, 



QCD sum rules 



28, 



29, 



381 ] involved in the SVZ sum rules. The basic idea of light-cone 



391 ] is to adopt the twist expansion of correlation functions near the light- 
cone instead of the dimension expansion of operators at short distance. Therefore, the essential inputs in 
the light-cone QCD sum rules is the hadronic distribution amplitudes other than vacuum condensates in 
the QCD sum rules. One important advantage of light-cone QCD sum rules is that it allows a systematic 
inclusion of both hard scattering effects and soft contributions [39J. In view of the above arguments, 



we will estimate the form factors for B c to charmonium states based on the light-cone QCD sum rules 
approach in this work. 

The structure of this paper is organized as follows: we first display the light-cone distribution am- 
plitudes of various charmonium states in section [Til The light-cone QCD sum rules for the form factors 
responsible for the decay modes B c — > M cS are derived in section IIIII The numerical computations of 
form factors in light-cone QCD sum rules are performed in section IIVI The decay rates for semileptonic 
decays of B c to various charmonium states, a brief analysis on comparisons with the results that obtained 
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with the help of other approaches in the literature and discussions on the S-D mixing of "0(3686) and 
-0(3770) in the weak decay of B c are also included in this section. The last section is devoted to our 
conclusion. 



II. THE LIGHT-CONE DISTRIBUTION AMPLITUDES OF CHARMONIUM STATES 



The light-cone distribution amp 
following non-local matrix element 



itudes (LCDAs) of pseudoscalar charmonium can be defined by the 
4C 



(P(p)\c(x) a c(0) p \0) = --Jp f 1 due iu v*[( l5 i>W» + mp(T6 (1) 



o 



4 

where <p v {u) and <p s {u) are twist-2 and twist-3 LCDAs of the pseudoscalar charmonium respectively. 
The decay constant fp can be determined generally by decay width of the double photons decay of the 



pseudoscalar meson as 



r(P^ 77 ) = 4 < 4 ™>^. (2 ) 



Making use of the branching fractions of r] c — > 77 and the full width of r] c 



BR(t/ c -> 77) = (2.8 ± 0.9) x 10" 4 , T Vc = (25.5 ± 3.4)MeV, (3) 

we can achieve the decay constant f Vc as 401^ygMeV. However, there is no data on 77" — > 77 till now, 
hence it is impossible to extract the decay constant of r]" directly from the experiments. In view of this 

43J 



point, we fix the decay constant f„" through the assumption, -^ £ - = 7^—, which has been used in 

/c JVC JJ/lp 

before. The decay constant of vector charmonium can be derived through leptonic decay V — > e + e~ as 



3m v T v ^ ee 

/y = V A^Ql ■ (4) 



42] 



Combining the above relation and the data given in 

Tj/^ ee = (5.55 ± 0.14 ± 0.02)keV, T^^ ee = (2.48 ± 0.06)keV, Y r ^ ee = (0.86 ± 0.07)keV, (5) 

we can obtain the decay constants of ijj(nS)(n = 1, 2, 3) as 

fj/^ = 416l^MeV, fa = 304llMeV, /^ = 187±8MeV. (6) 

In light of the assumption mentioned above, we arrive at the decay constant of rj" as 180^MeV. 
Moreover, the decay constant of rj' c can be determined as 293±|jjMeV. 

It needs to be pointed out that the tensor structure, which is suppressed in the heavy quark limit, has 
been neglected in the right hand side of the Eq. (pQ). When it comes to the explicit forms of 4> v (u) and 
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FIG. 1: The shape of the distribution amplitude 4> v (x) for 77" with v 2 = 0.3. 



4> s {u), we will adopt a simple model advocated in [45J]. Firstly, one should write down the Schrodinger 
equal-time wave function ^sch{f) for the Coulomb potential, and then perform the Fourier transformation 
of it to the momentum space as ^sch(k)- Next, in terms of the substitution assumption proposed in 46|] 
(see also Eq. (|A3P ). we can derive the expression of wave function ^ sch{ x ii from ^Sch(k)> where the 
momentum fractions xi, X2 of c and c quarks in the charmonium satisfy the relation x\ + X2 = 1. Finally, 
one can achieve at the LCDAs of charmonium $?sch(xi) by integrating over the transverse momentum 
k^. Based on this prescription, we can obtain the LCDAs for rj" as 



4> v {x) 



cj> s {x) 



10.8x(l 



x(l - x)[l - 4x(l - x)(l + fj)] 2 ^ l ~ v2 



2.1 



[l-4x(l - x)(l - 
x(l-x)[l-4x(l-x)(l + g)] 



(7) 

[l-4x(l-x)(l-f)P 

where the variable v reflects the mean charm quark velocity and is taken as v 2 = 0.30 ± 0.05 [45|] in the 
numerical analysis. To be more clear, the shape of the distribution amplitude <f> v (x) is shown in Fig. Q] 
with v 2 = 0.3. 

Similarly, we can also derive the LCDAs for rj' c 

x(l - x)(l - 2x) 2 



<^(x) 



<j>\x) 



10.6s 1 



[l-4x(l-x)(l-f)] 2 



2.1 



x(l - x)(l - 2x) 2 
[l-4x(l-x)(l-£)] : 



Similarly, the LCDAs of scalar charmonium state can be defined by 

(S(pMx) a c(0) p \0) = -Js f 1 due iu P- x U) 0aX v (u) +m s (I)/3aX s (u)}, 
4 Jo 



(8) 



(9) 



with x"( u ) an d X s ( u ) being the twist- 2 and twist-3 DAs for the scalar meson respectively. Based on the 
method of building the model for heavy quarkonium's distribution amplitudes given above, we can obtain 
the explicit forms of distribution amplitudes as 

x f x(l-x)(l-2x) 4 } l ~ v2 
X ix) = 90.2x(l - x (1 - 2x <^ - 2 — \ 

f x(l-x)(l-2x) 4 I 1 "" 2 . x 

X s {*) = 1-9 - \ - A Wh , Q , (10) 



[l-4x(l-x)(l-f)]3 



for the 2 Pq charmonium x' c o- l n addition, the decay constant of x' c o can be calculated as 2633 7 MeV 
making use of the assumption j-^- = and the value of f Xc0 , which was estimated to be 360MeV in 
Ref. 



4d, 



44]. 



The non-local matrix element associating with the vector charmonium can be decomposed as 

(V(p,e)\c(x) a c(0)p\0) = - f 1 due^ifymvij? - — fiV ± (u) + fymy— 6>V L {u) 

4 Jo I P ■ x p • x 



1 9m 1 

+$ 4fiV T (u) + 7 (fv -/y)W7 M 7 5 e*V^^(u) , (11) 

4 m V J 8a 



(3a 

where Vl(u), Vt(u) are the leading twist longitudinal and transverse LCDAs of vector charmonium, and 
V±( u )i Va(^) are the twist-3 ones. Following the methods described above, we can deduce the manifest 
expressions of these distribution amplitudes as 

s f x(l - x)((l - 2x) 2 [l - 4x(l - x)(l + £)} 2 "I l ~ v2 

V L {x) = 10.8x1 - x)\ -± 1K \ / { v2 ^ 

l [l-4x(l-x)(l-fg)] 4 J 

r / ^(x{l-x)((l-2x) 2 [l-4x{l-x)(l + l 4)} 2 ) 1 ~ v2 
V ± (x) = 1.74 1 + (1 - 2x) 2 }\ — — — j-^ ^IL\ 

I fl - 4x( 1 - x)(l - J 



F T (x) = Va(x) = V L (x), (12) 
for the 4 3 5'i charmonium. In the numerical calculations, the decay constants fy and fy are assumed to 



be equal [40|]. The purely leptonic decay of y(4260) — > e + e~ ie estimated to be 0.72keV 
we can obtain the decay constant /y(426o) as 176MeV. 

Similarly, the light-cone distribution amplitudes of 3 3 L>i charmonium can be derived as 

, ,f x 2 (l-x) 2 (l-2x) 6 ) l ~ v2 

V L (x) = 2.8x1 -x i — 4 — ^ 

l[l"4x(l-x)(l- J)P/ 

r N2nf X 2 (1-X) 2 (1-2X) 6 1 ^ 



211 ] , from which 



[l"4x(l-x)(l-g)P 
Vt(x) = Vk(x) = V L (x). (13) 



With the hypothesis ^| 3 ^ Pl) = y-^- and f^p^sD^ = 47.8 MeV [47(, we can achieve the value of f^,^^) 
as 21.5^Q5MeV under the above assumption. This is a quite small decay constant, comparing with that 



of the corresponding S-wave charmonium states. This will surely lead to the quite small form factors, 
since the transition form factors are proportion to the decay constant of the final state meson as can be 
observed form the light-cone sum rules in the next section. In the same way, we can arrive at the decay 
constant of ip(2 3 D 1 ) as f^ Dl ) = 34.9±{];fMeV. 

For the sake of investigating the 2S — ID mixing of "0(3686) and ^(3770), it is necessary to derive the 
light-cone distribute amplitudes for the tp^Si) and ^(l 3 -^) based on the model discussed above. To 
be more specific, the LCDAs for tp(2 3 Si) can be given by 

x(l-x)(l-2x) 2 ^- ! ' 2 



V L (x) = V T {x) = V A {x) = 10.6x(l - x) 



[l_4x(l-x)(l-^)]2 



4 

v ± ( X ) = i.7[i + (i - 2 X f)\ - x{ ]: x){ \: 2x)2 ^X ~ v2 ; (M) 

4 



[l-4x(l-x)(l-^ : 



while the LCDAs for ip(l 5 Di) can read as 

x 2 (l-x) 2 (l-2z) 2 ^- v2 



V L (x) = V T (x) = V A (x) = 3.6x(l - x) 



[l-4x(l-x)(l-^)]3 
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V x (x) = 0.77[1 + (1 - 2x) 2 ]{ ^ 2 ( 1 -^) 2 ( 1 -^) 2 Y~ V \ (15) 

l[l-4x(l-x)(l-^)] 3 J V ' 

Moreover, we also would like to present the explicit forms of LCDAs for ^(3 3 iSi) and tp(2 3 Di), which 
are essential to study the 3S — 2D mixing of "i/>(4040) and ^(4160). The LCDAs for ^(3 3 <Si) can be 
calculated as 

. ,'x(l-x)[l-4x(l-x)(l + ^1 
V L (x) = V T (x) = V A (x) = 10.8x(l - x)> — - - - 211 



^?\\ ^ ^ K \,: ; vJ )]2 Y~ v2 ; (i6) 



[l-4x(l-x)(l-f)P 

x(l -x)[l- 4x(l - x)(l + fj)] 2 
[l-4x(l-x)(l-^)P 

while the LCDAs for ip(2 3 Di) are given by 

V L (x) = V T (x) = V A (x) = 3.2x(l - x){ X J ~ X) \ %S ~ V2 

I 1 - Ax(l - x)(l - 4 J 



16 1 



V ± (x) = 0.70[1 + (1 - 2x) 2 }\ x2(1 X)2(1 2X) 1 V " . (17) 
V ' L V ' J \[l-4x(l-x)(l-fJ)] 4 J V ' 

As far as the axial- vector charmonium is concerned, the corresponding non-local matrix element can 
be analyzed as 



48] 



(A(p,€)|c(x) a c(0)^|0) = -- j 1 due iu v x \f A m A {,t*- — ^^(u) + f A m A — M\\(u) 

4tjQ ^ J) ' X J) ' X 



1 1m 1 

+fl 4 M±(u) + -AfA -fl)e^a^e*»p<*xPgl(u) , (18) 

4 m A J /3a 



where 4>\\(u), <fi±(u) are of twist-2, and g\(u) and gj_(u) are the twist-3 LCDAs of axial- vector charmo- 
nium. As for the states, (j>\\(u), g±(u) and gj_(u) are symmetric under the exchange of momentum 
fractions u and 1 — u, but </>±(u) is anti-symmetric under this exchange. On the contrary, <p±(u) is 
symmetric for n l P\ states, while 4>\\(u), g^(u) and gj_(u) are anti-symmetric in this case. Following the 
procedure of constructing the wave functions for heavy quarkonium shown above, we can arrive at 

1 1-v 2 

[1 - 4x(l -x)(l 



Mx) = 90.2x(l - - 2x){ g(1 ^ } 



for the 2 3 Pi charmonium, 



(j>Ax) = 10.6a(l-s){ X(1 2X)2 2 Y V , (20) 

V ' K ; l[l-4x(l-x)(l- ^)] 2 J V ' 

for the l 1 ^! charmonium h c , and 

, , x , xf x(l - x)(l - 2x) 4 l 1 "^ 2 , x 

= 9.5il-i i ^ — \ , 21 

y ' K ; l[l-4x(l-x)(l-^)]3j V ' 

for the 2 1 P\ charmonium h' c . 

As can be seen below, only the leading twist LCDA of axial- vector charmonium 4>±(u) is involved 
in the light-cone QCD sum rules of form factors, hence, the expressions of other three distribution 
amplitudes will not be shown here. Notice that we assume the decay constants Ja = /J in the practical 
numerical analysis, the same as that for the vector charmonium. Meanwhile, the decay constant of P wave 
charmonium 2 3 Pi(cc), was estimated to be 207MeV j^J very recently. Moreover, the decay constant of 
h r . and h' c are taken the same as that for the Xd m |3] and x'd respectively, namely fh c = f Xcl = 335MeV 



fK = = 207MeV. 



III. LIGHT-CONE QCD SUM RULES FOR THE WEAK TRANSITION FORM FACTORS 

For the semi-leptonic decays of B c — ► M cc -W[, the effective weak Hamiltonian is given by 

He f f(b clui) = ^V cb c^(l - y^bl-fQ- - 75 )^ + h.c. , (22) 

where V c b is the corresponding Cabbibo-Kobayashi-Maskawa (CKM) matrix element. In order to estimate 
the decay rates of B c — > M c5 Wi, we need to calculate the hadronic matrix element (M cg |c7 /i (l — ^)b\B c ) 
at first, which can be conventionally parameterized in the following forms: 

(Pcc(p)\5^b\B c (p + q)) = f + {q 2 )p^ + f-{q 2 )q^ (23) 
(5ce(p)|c7 M 766|S c (p + g)) = -i[U{q 2 )p„ + /-(<7 2 )<7/d\ (24) 



8 



(V cS (p)\c^(l - l5 )b\B c (p + q)) 



{A*(p)\&yJl-T5)b\B c (p + q)) 



m Bc + my 



e l iaf3~ ( e* a q (3 p' y - i(m Bc + my)A 1 (q 2 )e 



2\,* 



+ (e -g)(2p + ■ (e 



m Bc + my 
2iA{q 2 ) 

m Bc - m A 
V 2 (q 2 ) 



m Bc + my 



(25) 



2\,* 

A* 



(e*-<z)(2p + 



+ 



-(e* • ? )g MJ (26) 



where the anti-symmetric forth rank tensor is defined as Tr [7^7,^70-75] = 4ie /J ^p (J . 

Below, we will derive the general formulae for the form factors of B c — > M c5 (M = P, S, V, A) in the 



light-cone QCD sum rules approach. Following Ref. 



51 



521 ]. the correlation function is selected with the 



insertion of chiral current, to which the twist-3 distribution amplitudes of final states do not contribute 
at all for semi-leptonic B c — > P(S) decays. As for the B c — s> V{A) decays, the two-particle distribution 
amplitudes of twist-3 also have no effect on the correlation function with the insertion of chiral current 
in the heavy charm quark mass limit. Besides, the twist-3 distribution amplitudes relating to the three- 



particle ccg Fock state, which are suppressed by a factor (-^QCD/mcc) 2 45J] with m Ec being the mass 



of the charmonium, are also omitted in this work. The estimation of correlation functions in the QCD 
representation can be carried out following the standard prescription given in {53! . 5^ ]. 



A. Light-cone QCD sum rules for the weak transition form factors of B c — > P cE 

Based on the above analysis, we firstly construct the following correlator Tl fJ- (p,q) with the insertion 
of the chiral current: 

IL>,g) = i [ d i xe i ^{P c6 {p)\T{c{ X ) 1 Jl + 76)6(a;),6(0)i(l + 75)c(0)}|0). (27) 



One character of this correlation function is that twist-3 distribution amplitude of pseudoscalar char- 
monium has no influence on it and therefore the theoretical uncertainties can reduced considerably in 
this way. Inserting the complete sets of hadronic states with the quantum numbers the same as B c and 
making use of the following definition 

2 r 

(B c \U(l + 75 )c\0)= m y B ° , (28) 

we can arrive at the hadronic representation of correlation function (|27p as below: 

(Pcc(p)\c^(l + >y 5 )b\B c (p + q)){B c (p + q)\U{l + 7s)c|0) 



n M 0, q) 



m 2 Bc - (p + q) 2 
| y (P c5 (p)\c^(l+ l5 )b\h(p + q))(h(p + g)\bi(l+ 75 )c\0) 
k m 2 h -{p + q) 2 



(m b + m c ) (m 2 Bc - (p + q) 2 ) J s ^ c s - (p + q) 2 

where we have expressed the contributions from higher states of the B c channel in the form of dispersion 
integral with Sq c being the threshold parameter corresponding to the B c channel. On the other hand, 
we can also calculate the correlation function at the quark level: 

U,(p,q) = U^ CD (q 2 ,(p + q) 2 )p, + U QCD (q 2 ,(p + q) 2 )q, 

l Imng gD (s, g 2 ) [<*> l ImU QCD (s,q 2 ) 

ds I T \2 PV+ / ds / , ( 30 ) 

(m b +m c ) 2 VT S- (p + qY J(m b +m c ) 2 TT S - {p + q) Z 

Utilizing the quark-hadron duality assumption 

P Hs,q 2 ) = -lmUf CD (s, q 2 )e(s - s h ), (31) 

7T 

with i = "+, — " and performing the Borel transformation 

B M ,= lim tk±gj^ (')', (32) 



-( p+9 )2, n ^ 0O n! \d(p + qY 

-{p+q) 2 /n=M 2 

with variable (p + q) 2 to both two representations of the correlation function, we can finally derive the 
sum rules for the form factors 

2 m b + m c f s o c rrQCD, 2 n ( m B c ~ s \ , , QQ \ 

fi{Q) = —? — — / lmU i ( s 'g ) ex P — ur- ) ds - 33 

The QCD representation of correlation function (|27p can be calculated in terms of operator product 



expansion (OP 
transfer region 



54 



in both of the large space-like region (p + q) 2 <C — (mi, + m c ) 2 and the low momentum 



(nib — fTi c ) 2 — 2Agc , £)(mb — m c ) ~ 8.2GeV 2 , where the value of Aqcd is usually 
taken as 0.5 GeV. It is expected the light-cone QCD sum rules approach for the transition form factors 



will break down at large momentum transfer [5j], since the light-cone expansion for the description of 
final state meson is not well-pleasing in this case and the contributions from the higher twists would be 
important. The leading order contribution in the OPE can be gained simply by contracting the b-quark 



operators in the correlator (j27j) to a free b-quark propagator 



<0IH*)K0)|0> = / j^.—^, (34) 



which can be represented by Fig. [2] intuitively. 

Then we arrive at the correlation function at the quark level as 



f 1 4> v (u) 

U + (q 2 , (p + q) 2 ) = —2mbfp cS / du- — ^ [-contributions from higher twists, 

Jo (q + up) -ml + ie 

n_(g 2 , (p + q) 2 ) = + contributions from higher twists, (35) 

10 




FIG. 2: The tree level contribution to the correlation function Eq. (f2"T|) , where the current Js c (0) describe the B c 
channel and the current jp,(x) is associate with the b — > c transition. 



where the hig. 
charmonium 



ier twists contributions are at least from twist-4 distribution amplitudes of the pseudoscalar 



52, [55|, [5g] . Substituting the Eq. ([35l) to Eq. ([33j) . we can finally derive the light-cone QCD 
sum rules for the form factors fi(q 2 ) as below 



Mq 2 ) - 
Uq 2 ) = o, 



2m b (m b + m c )f Pct 



( m 



fB c m 2 B 



exp( 



B c 



M 2 



du 

— <p (ii)exp 
A u 



m b ~ u{q 2 — up 2 ) 
uM 2 



(36) 



up to the accuracy of twist-3 LCDAs, with 



A 



s — q 2 — p 2 ) + y (s — q 2 — p 2 ) 2 + Ap 2 (ml — q 2 ) 



(37) 



p 2 being the mass square for the corresponding charmonium state {m 2 P in this subsection) and s being 
the threshold value of B c channel. It needs to be emphasized that the vanishing of f~{q 2 ) up to the 
twist-3 LCDAs of pseudoscalar charmonium and leading order of the strong coupling constant a s is the 
consequence of the large-recoil symmetry [57]] , which emerges in the case of large recoil momentum for 



the final state meson and can be broken by the hard gluon corrections 



B. Light-cone QCD sum rules for the weak transition form factors of B c — > S c 



Following the derivation of the light-cone sum rules for B c — > P c5 , the correlation function of B c — > S C c 
can be written as 

H M (p,«) = i J d A xe ill x {S cS (p)\T{c(x) 7fl (l +Tfc)&(aO,&(0)i(l + 75 )c(0)}|0). (38) 

Matching the results of the above correlator calculated in the quark level and hadron representation 
respectively and performing Borel transformation with the variable (p + q) 2 , we can achieve the light- 
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cone sum rules for the transition form factors as below 



2m b {m b + m c )f Sc - c (m 



fs c m 2 



cxp 



M 2 



1 du v 

— X (^)exp 
A U 



m b ~~ u il 2 ~ U P 2 )~ 
~uM 2 



(39) 



where the lower limit of the integral A has been given in Eq. (|37D . 



C. Light-cone QCD sum rules for the weak transition form factors of B c — > V cS 

In the same way, the correlation function with the insertion of chiral current for B c — ► V C c can be 
chosen as 

H„(p,q) = i f d 4 x e ^(y c5 (p)|r{c(x) 7/ ,(l-75)fe(x),5(0)i(l+ 75 )c(0)}|0). (40) 



The hadronic representation of this correlator can be derived as 
u ( s _ (Vcc(p)\c7^ - l5)b\B c (p + g))(B c (p + q)\U(l + 75 )c|0) 

1 W<1) - _2 {p + q) 2 



+E 



B c 

(Vcc(p)\c^(l - l 5 )b\h(p + q))(h(p + q)\bi(l + 75 )c|Q) 
m h ~ (P + ?) 2 

^m 2 B J Bc V(q 2 ) w p 7 _ m 2 B J Bc (m Bc + m Ve ,)Ai(g 2 ) # 



(m b + m c ){m Bc + m Vc -)(m 2 Bc - (p + g) 2 ) 7 (m 6 + m c ){m 2 Bc - (p + g) 2 ) 

| ; m| e / Be (e* • g)l4 2 (g 2 )(2p + g)„ + ^(g 2 )^] + j°° ^ p h v (s,q 2 ) ^ ^ ^ 
{m h + m c )(m Bc +m Vc c){m 2 Bc - (p + q) 2 ) J s Bc s - (p + q) 2 A 



DO 



Jfr + J.? s-(p + q) 2 [e Q} - (4ij 

Besides, the correlation function in Eq. (l40p can also be formulated as 

IL>, q) = n® CD (q 2 , (p + g) 2 )6^ 7 e w 9 V - »H§f V, (p + g) 2 )< 

+iII« C V, (p + g) 2 )(e* • g)(2p + g)„ + iU^ D (q 2 , (p + g) 2 )(e* ■ q)q, 



(mi,+m c )2 VT S - (p + g) 2 J(m b +m c ) 2 VT S - (p + 



2 A* 



oo 



1 Imn«f( S , g 2 )(2p + g) M + Imn« CD ( S) g 2 )g M . 



J{m b +m c y vr s - (p + q) A 

Matching these two representations of the correlator and performing the Borel transforming with the 
variable (p + g) 2 on them, we can obtain the light-cone QCD sum rules for the form factors of B c — > V C c 
as 



(m b + m c )(m fl +m v j lmJl QCn 2 frr^ 

2nf Bc m% J(m b +m c )* v \ M 2 
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171b + m c 



nfB c m%(m Bc + m Vc - c ) J(m b +m c ) 



lmUf D (s,q 2 )e W (^ 7 ^)ds, 



(m b + m c )(m Bc + my cE ) f s o 



nfB c m B 



Im Il ( ^ D {s,q z )eyiY> 



(m b +m c ) 2 



M 2 

m B — S 
M 2 



{m b + m c )(m Bc + m Vc5 ) f s o 



Imn A3 (s,q 2 )expi 

(rnb+rric)^ \ J-Vl 



m 



B<- 



ds, 



ds. 



(43) 



^fB c m 2 Bc 

Substituting the QCD representation of the correlation function in Eq.(|40p with the help of the OPE 
technique, we can derive the explicit forms of the form factors in the light-cone QCD sum rules as 

m b ~ u (l 2 — up 2 ) 1 



V(q 2 
Mq 2 



V M A ) J a u 



fB c m 2 



u 



M 2 



(m b + m c )f VcE f m B c \ f 1 du 

JB c m Bc (m Bc + m VcB ) \M Z J J A u 
A 2 (q 2 ) = -A 3 (q 2 ) = V(q 2 ), 



m 2 — u(q 2 — up 2 ) 1 m 2 — q 2 + u 2 p 2 



iM 2 



(44) 



with the lower integral limit A defined by the Eg. (1371 ). It needs to note that similar results were also 
obtained in Ref. [561 . 159] . 



D. Light-cone QCD sum rules for the weak transition form factors of B c — > A c 



The derivation of light-cone QCD sum rules for B c — > A c5 is very similar to that for B c — ► V C c discussed 
before. The correlator for B c — > A c5 can be given by 

n„(p, q ) =i J d 4 xe^(^ cc -(p)|r{c(x)7 / ,(l-75)^),K0)^(l + 75)c(0)}|0). (45) 

We will skip the detailed derivation of sum rules for the form factors in B c — > A c5 and only display the 
final results of them as 



A(q 2 ) 

Vi(q 2 ) 
V 2 (q 2 ) 



(m b + m c )(m Bc - m AcE )f Ac 



fB c m 2 
(m b + m c )fj 



B, 



B, 



da 



m b — u {q 2 — u p 



2\i 



U 



M 2 



f Bc m 2 Bc (m Bc - m Ac5 ) 
-V 3 (q 2 )=A(q 2 ). 



cxp 



m, 



M 2 



du 



A U 



4>±(u)exp 



ml — u{q 2 — up 2 ) 1 



uM 2 



2^2 



mi — q + up 



(46) 



13 



IV. NUMERICAL RESULTS FOR THE FORM FACTORS AND DECAY RATES 



Now we are going to analyze the sum rules for the form factors numerically. Firstly, we collect the 



input parameters used in this paper as below 



42. 



61 



62 



63 



64J]: 



m e = 


0.511MeV, 


m T 


= 1.777GeV, 


m b = 


(4.68 ± 0.03)GeV, 


m c 


= (1.275 ±0.015)GeV 


m Bc -- 


= (6.286 ±0.005)GeV, 


Ib c 


= (395 ± 15)MeV, 


T B C = 


(o.463±g:8g)i», 


b c 
V 


= (45 ± l)GeV 2 , 


Gf = 


1.166 x 10- 5 GeV" 2 , 


\v cb 


= (42.21±8:S°) x io- 



(47) 



It is noted that the decay constants of various charmonium states have been discussed comprehensively 
in section [TT1 

The choice of the threshold parameter s can be determined by the condition that the sum rules should 
take on the best stability in the allowed M 2 region. Besides, the value of threshold parameter should be 
around the mass square of the corresponding first excited state, hence they are also chosen the same as 
that in the usual two-point QCD sum rules. The standard value of the threshold in the X channel is 
sqx = (mj + Ax) 2 , where Ax is usually taken as 0.5GeV 



65, 



66 



67, 



681 ] approximately in the literature. 



To be more specific, we will adopt the threshold parameter forl? c channel Sq c as (45 ± l)GeV 2 for the 
error estimate in the numerical analysis as shown above. 

It is well known that the form factors should not depend on the Borel mass M in the complete theory. 
However, we can only truncate the operator product expansion up to some finite dimension and perform 
the perturbative series in a s to some order in practice, both of which will result in the dependence of the 
form factors on the Borel parameter definitely. Therefore, one should find a region where the results only 
depend moderately on the Borel mass, and the approximations for the above truncations in the complete 
theory are reasonable and acceptable. 

In general, the Borel mass M should be chosen under the requirement that both the contributions 
from the higher resonance states and higher twist distribution amplitudes are small (no more than 30 
%) to ensure the validity of the OPE near the light-cone and the quark-hadron duality being a good 
approximation. As for the decay of B c — ► rj' c , we indeed find the Borel platform M 2 £ [20,30]GeV 2 , 
which is also consistent with the number obtained in the two-point QCD sum rules corresponding to 
the decay constant of fg c 69J. The light-cone QCD sum rules of form factor f+{q 2 ) at zero momentum 
transfer is shown in Fig. [3j The values of /+(0) with various uncertainties rooting in Borel mass, threshold 
value, decay constants of the related mesons, heavy quark masses and the parameter v 2 involved in the 
LCDAs of charmonium have been collected in Tabled! from which we can find that the total uncertainties 
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TABLE I: The form factors fi(0) responsible for B c — > P(S) c5 decay in the light-cone QCD sum rules approach; 
the errors for these entries correspond to the uncertainties in the Borel mass, threshold value, quark masses, decay 
constants of two mesons and variations of v 2 in the LCDAs of charmonium respectively. 



decay modes 


/+(0) 


/-(0) 


B c - ^(2* So) 


n 09+0.03+0.02+0.01+0. 17+0.01 

u - oz -0. 01-0. 02-0. 01-0. 19-0. 01 





B c -> X(3940)(3 1 S , ) 


n AR+°- 01+0. 00+0. 00+0. 10+0. 01 
u - w -0.01-0.01 -0.01-0. 11-0.01 





B c F(3940)(2 3 P ) 


O c+0. 1+0. 0+0. 0+0. 2+0.0 
z -°-0.1-0.1-0.1-0.2-0.2 






of form factors are indeed at the level of (20 — 30)% as expected by the general understanding of the 
theoretical framework. The form factor /-(0) up to the twist-3 LCDAs of P c g and leading order of a s 
is zero as a result of the large-recoil symmetry. The q 2 dependence of the form factor f + (q 2 ) calculated 
from light cone sum rules is shown in Figf4]in the physical kinematical region < q 2 < (rriB c — m^) 2 - 



Since the number of (m# c — m^/J 2 ~ 7.0GeV 2 with m^i = 3.638 ± 0.004GeV 42] being used, is smaller 
than that of (m& — m c ) 2 — 2AQCD{ r nb — m c) — 8.2GeV 2 , the OPE technique near the light-cone can be 
performed in the whole kinematical region effectively. 

We can further evaluate the sum rules for the form factors associating with B c to other charmonium 
states. For example, the only difference for the calculation of decay mode B c — ► X(3940)(3 1 S'o)/^, is 
to substitute the LCDAs of X(3940)(3 1 S'o) for that corresponding to r}' c compared with the decay of 
B c — > if c Wi. In the light of Eq. (|39p and the light-cone distribution amplitudes of scalar charmonium 
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calculated before, it's straightforward to estimate the light-cone sum rules for the transition form factors 
of B c — > x'co^lf the number of which has been grouped in Table HI Evaluations of the form factors relating 
to the B c — ► V{A){cc) decay are also easily carried out with the help of Eq. (|44|46p and the LCDAs of 
(axial) vector meson displayed in Section [TIJ Since the calculations for all of these form factors are quite 
similar, we will not explicitly repeat the details anymore and only display the final results in Table HT1 

Utilizing the above form factors and the input parameters shown in Eq. (|47p . we can proceed to 
compute the branching ratios of these modes. Followin g th e standard procedure, the differential partial 
decay rate for B c — ► M C clvi (I = e, t) can be written as [42j 



dg 2 " (2tt)3 32m| c J Umi 



/"Umax — _ 

/ \M Bc ^ McdPl \ 2 du, (48) 

«/ Urn.i-n. 



where u = (pm cE +w) 2 an d Q 2 = (pi +Pp ( ) 2 ; Pm^i Pi an d pp t are the momenta of M c5 , I and v\ respectively; 
M is the decay amplitude after integrating over the angle between the I and M c5 . The upper and lower 
limit of u are given by 



u m ax = (E* McB + Eff - (^E* A l_ c - m 2 Mce - ^Ef - m 2 ) 2 , 

Umin = (E* Mce + Et) 2 - {^Ef Ic _ c -m 2 Mc _ c + ^Ef-m 2 ) 2 ; (49) 

where E^ _ and E£ are the energies of the charmonium state and the lepton in the rest frame of lepton- 
neutrino pair respectively and the manifest expressions of them can be given by 



m B c ~ m ll c - c ~ <f p * _ q 2 +m 2 

m C c irj^ > - n m • l 50 J 

q 
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TABLE II: The form factors /i(0) responsible for B c — ► V c5 (A c5 ) decays in the light-cone QCD sum rules approach; 
the errors for these entries correspond to the uncertainties in the Borel mass, threshold value, quark masses, decay 
constants of two mesons and variations of v 2 in the LCDAs of charmonium respectively. 



decay mode 


7(0) 


A 1 (0) 


B c -> </>(2 3 Si) 


n Qn+o° 3 +o-02+o. 02+0. 04+0. 01 
u - JU -o. 02-0. 03-0. 02-0. 05-0. 00 


n 00+0.01+0.01+0.00+0.02+0.00 

u - J °-0. 01-0.02-0.01-0.02-0.00 


B c ^Dt) 


n 11+0.00+0.01+0.00+0.00+0.00 
u - 11 -0.00-0. 01-0. 00-0. 00-0. 01 


A q+0. 0+0. 5+0. 2+0. 2+0.1 1n -2 
^• y -0. 1-0.6-0.3-0.2-0.3 x iu 


B c -> ?/>(3 3 Si) 


n r9+0. 02+0. 00+0. 00+0. 04+0. 00 
u - oz -0. 01 -0.01-0. 01-0. 04-0. 00 


n 91+0.01+0.00+0.00+0.02+0.00 

U - ZL -0.01 -0.00-0.00-0.02-0.00 


B c -> V>(2 3 £>i) 


7 9+0.0+0.9+0.5+0.4+0.6 v 1 n -2 
' - z -0. 1-1. 0-0. 6-0. 5-0. 6 x iu 


n +0. 0+0. 4+0. 2+0. 2+0. 3 v 1 n -2 
J ' u -0.0-0.4-0.2-0.2-0.2 A iu 


B c -► r(4260)(4 3 S' 1 ) 


n 47+0.01+0.00+0.00+0.07+0.00 

u '^' -0.01-0.01-0.01-0.07-0.01 


n 1 0+0. 01+0. 01+0. 00+0. 03+0. 01 

u - io -0. 00-0. 00-0. 00-0. 03-0.0 


B c ^V(3 3 -Di) 


0+0. 0+0. 6+0. 3+0. 3+0. 4 1 n -2 
°'°-0.1— 0.7-0.4— 0.3— 0.5 A iu 


1 r+0.1+0.3+0.2+0.1+0.2 ln -2 
L -°-0. 0-0. 2-0. 1-0.1-0.1 A iu 




A(0) 


K(0) 


P c -> X(3872)(2 3 Pi) 


n r;o+0. 02+0. 01+0. 01+0. 02+0. 04 
u - oo -0. 03-0. 00-0. 00-0. 02-0. 03 


7R +0. 12+0. 05+0. 02+0. 14+0. 24 
J - ' D -0. 18-0. 02-0. 00-0. 14-0. 20 


P c -> X(3940)(2 3 Pi) 


n n +0.02+0.01+0.01+0.02+0.04 
u - di -0.02-0. 00-0. 01 -0.02-0. 03 


07+0. 11+0. 05+0. 02+0. 15+0. 24 
-0.19-0.03-0.02-0.15-0.22 


^ -> ^(^Pl) 


n 9 o+0.01+0.01+0.01+0.01+0.00 
u - zo -o. 01 -0.00-0. 00-0. 01-0. 00 


1 r i +0.04+0.04+0.02+0.06+0.01 
L - 01 -0.02-0. 04-0. 02-0. 06-0. 00 


P c -> ^(2^1) 


n -14+0. 00+0. 00+0. 00+0. 01+0. 00 
u - iq: -o. 00-0. 01-0. 00-0. 01-0. 00 


1 i n +0. 00+0. 00+0. 04+0. 00+0. 00 
1 - iu -0. 00-0. 00-0. 04-0. 00-0. 00 



The numerical results are shown in Table IIIII and IIVI together with the numbers obtained in other 
approaches for comparison. It is observed that the decay rates for B r — > hr and Br — > r/' c calculated 



61 



m 



71 



72 



731 ] within the error 



in this work are consistent with that obtained in other frameworks 
bars, such as quark model, Bethe-Salpeter equation, SVZ sum rules and so on, therefore, the branching 
fractions for B c to the new charmonium states presented in this work are reliable and acceptable. 

It should be noted that the assignment of X(3940) as 3 1 Sq charmonium state leads to the production 
rate as 1.9 x 10 -4 in the weak decay of B c — ► X(3940)eP e , the magnitude of which is one order smaller 
than that for the interpretation of X(3940) being a 2 3 Pi charmonium. This particular phenomenology 
can provide valuable information for us to discover the inner structures of X(3940). Besides, both of 
X(3872) and Y(3940) should be observed in the weak decay of B c -> X(3872)/Y(3940)e^ e in view of 
the branching ratio as large as 10~ 3 order, on the condition that they can be explained as 2 3 Pi and 
2 3 Pn charmonium states respectively. If future experimental measurements deviate from our predictions 
heavily, it will rule out the current J PC assignments of the charmonium states. It also needs to mention 
that the decay rates for semideptonic decays of B c — > M C c,tv t are also displayed in Table ITVl from which 
we can find that they are about one order smaller than the corresponding channel B c — > M C cev e due to 
suppression from the phrase space and sensitive dependence of form factors on the momentum transfer 
q 2 . In particular, the branching fraction of B c — > Y(4260)tP t is about two orders smaller than that for 
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TABLE III: Branching fractions of B c — > M c5 eD e semi-leptonic decays in the light-cone QCD sum rules approach; 
the errors for these entries correspond to the uncertainties in the Borel mass, threshold value, quark masses, decay 
constants of two mesons, lifetime of B c and variations of v 2 in the LCDAs of charmonium respectively. 



decay modes 


BR(this work) 


Other works 


B c -> V ' c (2 1 S )e9 e 


i 1+0.1+0.0+0.0+0.4+0.4+0.0 v in -3 
1 - i -0. 0-0. 0-0. 0-0. 5-0. 4-0.0 * iu 


3.2 x 10~ 4 [70] 
5.1 x 10" 4 [71] 


B c -> X(394O)(3 1 S )eu e 


i q+0. 2+0. 1+0. 0+0. 8+0. 7+0.0 i n -4 
J-^-0. 1-0. 1-0. 0-0. 9-0. 7-0.0 * lu 




B c -> F(3940)(2 3 P )ez> e 


7 9+0.9+0.1+0.1+0.9+2.5+0.0 i n -3 
' ' z -0. 5-0. 0-0. 0-0. 9-2. 8-0. 9 * iu 




B c -> y(4260)(4 3 S , i)eP e 


i r+0. 1+0. 0+0. 1+0. 1+0. 5+0.0 1n _4 
i,o -0. 1-0.1-0. 0-0. 1-0. 6-0.0 x lu 




P c A(3872)(2 3 Pi)eP e 


f- 7+0.9+0.0+0.1+0.5+2.3+0.7 -. n -3 
°-' -0.5-0.0-0.0-0.5-2.6-0.7 A lu 




P c -> A(3940)(2 3 Pi)ei? e 


fi n +0. 7+0. 0+0. 0+0. 5+2. 1+0. 6 v i n-3 
°- u -0. 5-0. 1-0. 1-0. 5-2. 3-0. 7 A lu 




P c h c {l x P l )ev e 


9 n+0. 3+0. 0+0. 0+0. 2+1. 0+0.1 v i n -3 
z ' J -0. 1-0.0-0.0-0.2-1. 1-0.0 A lu 


2.7 x 10" 3 [61] 
1.7±g;2 x 10~ 3 [72] 


B c -» h' c (2 1 P 1 )eP e 


K o+0. 3+0. 1+0. 0+0.4+1. 8+0.1 in -4 
d -°-0.2-0. 1-0.0-0.4-2. 1-0.0 A lu 





TABLE IV: Branching fractions of B c — > M cS tD t scmi-lcptonic decays in the light-cone QCD sum rules approach; 
the errors for these entries correspond to the uncertainties in the Borel mass, threshold value, quark masses, decay 
constants of two mesons, lifetime of B c and variations of v 2 in the LCDAs of charmonium respectively . 



decay modes 


BR(this work) 


Other works 


B c i^S^TVr 


o 1+0.9+0.1+0.1+3.3+2.8+0.1 v i n-5 
°- i -0. 5-0. 1-0. 1-3. 7-3. 2-0.0 x iu 


1.6 x 10" 5 [73] 


B c -> AT(3940)(3 1 S , )ri> T 


r 7+0.6+0.7+0.3+2.4+2.0+0.0 v i n -6 
d - ' -0.3-0.4-0.3-2.7-2.2-0.1 * iu 




B c -v y(3940)(2 3 P )-rz? T 


9 7+0.4+0.0+0.0+0.3+0.9+0.0 v i n-4 
z - ' -0.2-0.0-0.0-0.3-1.1-0.3 A lu 




B c -> r(4260)(4 3 S' 1 )ri7 T 


f. 1+0.5+0.8+0.3+0.5+2.2+0.1 1n -7 
°-^-0. 3-0. 4-0. 2-0. 5-2. 5-0.0 x iu 




B c -> AT(3872)(2 3 Pi)tz? t 


o 9+0.5+0.0+0.0+0.2+1.1+0.4 , n -4 
°- z -0. 2-0. 2-0. 0-0. 2-1. 3-0. 3 x iu 




P c -> AT(3940)(2 3 Pi)t^ t 


9 9+0.3+0.0+0.1+0.2+0.8+0.2 i n-4 
z - z -0. 2-0. 0-0. 0-0. 2-0. 9-0. 3 x lu 




P c -> ^(l^jTMr 


o 7+0.4+0.1+0.1+0.3+1.3+0.1 i n-4 
°- ' -0.2-0.1-0.0-0.3-1.4-0.0 A iu 


1.7 x 10" 4 [61] 
l-5±8:o >< W- 4 [72] 


B c ^h' c {2 l P 1 )rD T 


9 n +0. 2+0. 0+0. 0+0. 2+0. 7+0.1 , n -5 
z - u -0. 1-0. 0-0. 0-0. 2-0. 8-0.0 x iu 





the B c — > y(4260)ez^ e mode, since the sum of the mass for y(4260) and r lepton is almost close to the 
threshold of B c meson. 

Finally, we are in a position of concentrating on the S — D mixing of various vector charmonium. 
It's known that the S — D mixing of ^(3686) and -0(3770) may be essential to explain the large leptonic 
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decay width of "0(3770), the notorious pir puzzle 



47J and the enhancement of -0(3686) -> K L K S [7J]. The 



production of -0(3770) in B meson decays B + — > 0(377O)i^ + is found to be surprisingly large by Belle 
75 1 . which can be even comparable to B + — ► ip(3686)K + \t\. Hence, it is helpful to investigate the 
weak production of 0(3686) and -0(3770) in B c decays in order to test the above mixing scheme further 
and clarify the inner structures of them. Assuming that the physical state -0(3686) and 0(3770) are the 
mixture of l 3 _Di and 2 3 5i states, we have 

|0(3686)) = cosfl|2 3 5i) + sin#|l 3 £i), 

|0(377O)) = -sin#|2 3 Si) + cos#|l 3 L>i). (51) 

As for the mixing angle 9, two solutions 9 = — (12 ± 2)° or 9 = (27±2)° Q, Q]' were found in order to 
reproduce the leptonic widths of -0(3686) and -0(3770) [471 ] . The small mixing solution, i.e., 9 = — (12±2)°, 



is consistent with couple-channel estimates 
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8l| and the El transition -0' — > XcJl 
Based on the transition form factors of B c — > -0(2 3 5i) and B c — ► ^(l 3 Z?i) listed in Table [Til we 
can plot the production rates of them in the B c decays as functions of the mixing angle 9, which are 
displayed in Fig. [5] As for the favored mixing angle 9 = —12°, the branching fraction of B c — > tfj (3686) ei7 e 
is 1.5 x 10~ 3 , which is almost the same as the number of 1.7 x 10 -3 in the case of 2 3 5i state without 
mixing. However, the mixing component of i/j(2 3 Si) in the structure of ^(3770) is in particular important, 
which can increase the decay rate of B c -> -0(377O)ez7 e from 4.5 x 10" 5 to 2.1 x 10~ 4 . The reason is that 
the decay constant of 1 3 D± charmonium (47.8 MeV) is too small compared with that of 2 3 Si state (304 
MeV), therefore, even a small mixing angle can affect the decay rate of B c — ► ^(3770)ei^ e drastically, but 
almost has no effect on the decay B c — > -0(3686)ez7 e . 

Besides, exploring the properties of Y(4260) and ^(4415) as the mixing of 45 and 3D states in the 
weak decays can also shed light on the the inner structures of these charmonium-like particles. On the 
one hand, the assignment of 3 3 -Di charmonium for ^(4415) is supported by the fact that -0(4415) is 
dominated by the decay of ^(4415) — > DD^^AQQ) reported by the Belle Collaboration very recently 



82J. On the other hand, Y(4260) can be accommodated as 45 state naturally based on the analysis of 



211 ]. Moreover, the absence of Y(4260) signal in e + e — ► hadrons 

. Similar to the 5 — D mixing of ^(3686) 



production and decay characters of it 

can be explained quantitatively in the 5 — D mixing scheme 
and ^(3770) , we express the states of Y(4260) and ^(4415) as 

| Y (4260)) = cos#|4 3 5i) + sin0|3 3 L>i), 

|^(4415)) = -sin0|4 3 5i) + cos9\3 3 D 1 ). (52) 

In the above mixing picture, we can analyze the dependence of production rates for Y(4260) and 
-0(4415) in the weak B c decays on the mixing angle 9 with the help of the form factors of B c — > ^(4 3 5i) 
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FIG. 5: Decay rates of B c — > ^>(3686)ez7 e and B c — > ^(3770)ei7 e as functions of the mixing angle 9. The solid line 
represents the case of B c — ► ?/'(3686)ez7 e , while the dashed line is for B c — ► ?/>(3770)ez7 e . 




represents the case of B c — > F(4260)e^ e , while the dashed line is for B c — ► -0(4415)e^ e . 



and -B c — > "i/;(3 3 -Di) calculated before. As shown in Fig. [6l it can be observed that the decay rate for 
B c — > y(4260)eP e and B c — > ip(4A15)eu e are 1.5 x 10~ 4 and 1.1 x 10~ 6 for the null mixing angle. As a 
simple test, we find that the production rate of ^(4415) in B c decay can reach as large as 1.0 x 10~ 5 for 
the mixing angle of 6 = —12°; while the branching fraction of B c — > Y(4260)eP e is 1.4 x 10 -4 , almost the 
same as that in the case with zero mixing angle. 

For the completeness, we also consider the ?/>(4040) and ^(4160) being the mixing states of 3S and 
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FIG. 7: Decay rates of £? c — > ^(4040)ez7 e and £? c — > ^(4160)ei7 e as functions of the mixing angle 9. The solid line 
represents the case of B c — ► ^>(4040)ez7 e , while the dashed line is for B c — ► T/>(4160)ez7 e . 



The dependence of production rates for B c — > -0(4040)6^ and B c — > -0(416O)ez7 e on the mixing angle 6 
have been plotted explicitly in Fig. Without the mixing of S — D states, the branching ratios for 
B c — ► V ; (4040)e;? e and B c — > ■0(416O)e^ e are 2.9 x 10 -4 and 7.1 x 10~ 6 respectively. As for the mixing angle 
6» = -12°, the decay rates become 2.6 x 1CT 4 and 3.3 x 10~ 5 for B c -> ^(4040)e^ e and 5 C -> V(4160)ez? e , 
from which we can observe a large enhancement for the production of ^(4160) as we expect. 

Up to now, we only focus on the discussions of S—D mixing for the decays of B c — > V C cev e , which can be 
readily generalized to the case for the B c — > V C ci~v T decays. Subsequently, similar observations to the final 
states being the e v e pair can be achieved: Even a small mixing angle can result in considerable effectrs 
on the decay rates of B c — ► i\)(rc , D\)Tv T (n = 1,2,3), while branching fractions of the corresponding 
channels B c — > tp((n + 1) 3 <Si)tP t do not vary significantly. 

V. SUMMARY 

A number of new heavy charmonium states, such as rj' c , h c , X(3940), Y(3940), X(3872) and y(4260) 
are observed during the past several years. There exist various explanations for the quark components 



2D as 



1-0(4160)) 



1-0(4040)) 



cos0|3 3 Si) +sin#|2 3 L>i), 
-sin#|3 3 Si) + cos#|2 3 L>i). 



(53) 
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of the heavy mesons Y(3940), Y(3940), Y(3872), Y(4260) so far, such as charmonium states, tetraquark 
pictures, molecular bound states and so on. It is still early to give a definite answer for their solutions. 

In this work, we mainly focus on the charmonium interpretation of all the states, rj' c , h c , h' c , Y(3940), 
Y(3940), Y(3872), and Y(4260) produced in the exclusive semi-leptonic weak decay of B c meson. In 
order to compute the branching ratios of semi-leptonic weak decays of B c , we need to deal with the 
hadronic transition matrix element {M c5 \j^\B c ) , which defines the form factors governed mainly by non- 
perturbative QCD effects. In this paper, the light-cone QCD sum rules approach is used to evaluate 
various form factors. We choose the correlation function with the insertion of chiral current following the 



Ref. 
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52l | . the consequences of which are that the twist-3 LCDAs do not contribute to the sum rules 



for B c decays to the pseudoscalar charmonium and also for B c decays to the (axial) vector charmonium 
in the absence of three-particle wave functions. 

With the help of form factors calculated in the light-cone QCD sum rules approach, we give the decay 
rates for semi-leptonic decays of B c — > h c and rj' c , which agree with that derived in other frameworks. 
Besides, it is found that different interpretations of X(3940) can result in remarkable difference of the 
production rate in the B c decays, which would help to clarify the quark structures of the X(3940) with 
the forthcoming LHC-b experiments. Furthermore, the weak productions of X(3872) and Y(3940) in 
B c decays are large enough to be detected in the future experiments, supposing that they are indeed 
2 3 Pi and 2 3 Po charmonium states respectively. It is also observed that the mixing component of 2 3 S\ 
charmonium state in the structure of ^(3770) can enhance its production rate in B c decays heavily, even 
for a small mixing angle. Besides, the production character of Y(4260) and ^(4415) as the mixing of 45 
and 31? states as well as V(4040) and ^(4160) being the mixing states of 35" and 2D are also included 
in this work. In fact, all these decay rates depend heavily on the J PC assignments of the charmonium 
states. Therefore, our calculations can be used in LHC-b experiment to explore the components of these 
hidden charm mesons. 
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APPENDIX A: AN EXAMPLE OF CONSTRUCTING LCDAS FOR C HARMONIUM STATES 



re construction of LCDAs for heavy 



40 



45| described in section HH 



Taking rj" meson as an example, we would like to explain t 
quarkonium step by step in this appendix based on the procedure 

Firstly, we write down the radial Schrodinger wavefunction of n = 3, I = state for the Coulomb 
potential as 



ipSch(r) oc [1 - ^q B r + J^(gsr) 2 ]exp(-^), 



(Al) 



where qs is the Bohr momentum. Performing the Fourier transformation of the above wavefunction, 
then we can arrive at 

q% - 30q 2 B k 2 + 81/c 4 



(A2) 



(9k 2 + q 2 B ) 4 ' 

with k 2 being the square of three momentum, namely k 2 = |k| 2 . In terms of the substitution assumption 
46] 



(2x - 1) 



m 



m 2 + k 2 ; 



rrin 



2 ' " u x(l-x) 
we should make the following replacement towards the variable k 2 

, 2 k 2 ± + (l-2x) 2 m 2 



(A3) 



k 



4x(l - x) 

Now, we can derive the Schrodinger wavefunction for rj" as 

^Sch(x) oc J d 2 k_L^ 5ch (x,k ± ) 

J x(l-x)[l-4x(l-x)(l + ^)] 2 

OC X(l — X)< 5 — 

I [l-4x(l-x)(l-ftp 
where v 2 is defined as v 2 = q^/m 2 . Following the Ref. Q, we propose the LCDAs of rj" as 



(A4) 



(A5) 



^ a {x) = ^ y {x) 



x(l -x)[l- 4x(l - x)(l + §^)] 2 ^ 



[l-4x(l-x)(l-^)P 



(A6) 



where the power 1 — v 2 reflects the small relativistic corrections to the Coulomb wavefunctions 45(. It 
can be observed that above distribution amplitudes have the correct asymptotic behavior for both the 
heavy quarkonium in the heavy quark limit v 2 — ► and light mesons in the v 2 — ► 1 limit. 
Moreover, it is known that the asymptotic forms of pseudoscalar mesons can be given by 



4> V asy(x) CCX(1- X), C„W « 1. 



asy v 



(A7) 
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Substituting Eq. (|A7j) into Eq. (|A5|) . we can obtain the LCDAs for 77" 

' x(l -x)[l- 4x(l - x)(l + fij)] 2 ^ ^ 



0"(a:) = 10.8x(l 



[1 - 4x(l - s)(l 



_ 21 r I ( 1 - l)[1 -4,( 1 _ (A8) 

I [1 - 4x(l - x)(l - ^)] 3 J 

corresponding to the given value of v 2 = 0.30 as displayed in the text, where the normalization condition 
Jq cj) v ' s (x)dx = 1 has been used in the derivation of above LCDAs. In addition, we find that the fluctu- 
ations of the phenomenological parameter v 2 do not have significant effects on the shape of distribution 
amplitudes for charmonium states generally within the acceptable range v 2 = 0.30 ±0.05, which can also 
be verified from the numbers of decay rate for B c — > M cS Wi grouped in Table IIIII and IIV1 
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